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Abstract 

We present results of our simulations of QCD with two light dynamical quarks on a 32^ x 8 lattice 
at a current quark mass tuned to have the Goldstone pion mass of about 230 MeV. Employing 
the Taylor expansion method we proposed earlier, we estimate the radius of convergence of the 
series for the baryonic susceptibility by using terms up to eighth order Together with our earlier 
results, corresponding to the same physical parameters but on coarser lattices at respectively 1.33 
times and twice the lattice cut-off (a), we were able to attempt a march towards the continuum 
limit. 
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1. Introduction 

Whether the phase diagram of strongly interacting matter, governed by Quantum Chromo 
Dynamics (QCD), has a critical point in the temperature {T) and baryon chemical potential 
plane, is an exciting question that has attracted many theorists, phenomenologists, and heavy 
ion experimentalists. A variety of models, which have been successfully tested for hadronic 
properties in our world, such as an effective chiral Nambu-Jana Lasinio type model, lead to a 
phase diagram 1 1 ] with a critical point in a world with two light quarks and one heavier quark. It 
is clearly desirable to obtain it from QCD directly, or show a lack of its existence, especially since 
enormous efforts at RHIC and other accelerators are being devoted to look for it. One usually 
has to deal with large coupling constants in the world of (low energy) hadronic interactions. 
Non-perturbative lattice QCD, defined on a discrete space-time lattice, has proved itself to be 
the most reliable technique for extracting such information from QCD. The hadron spectrum 
has been computed successfully and predictions of weak decay constants of heavy mesons have 
been made. Since these fix quark masses and Kgco, a completely parameter free application of 
this approach to finite temperature QCD has yielded a slew of thermodynamics determinations, 
such as the pressure as a function of temperature. It is therefore natural to ask whether lattice 
QCD can help us in locating the QCD critical point. 
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Due to the well-known fermion doubling problem, one has to make a compromise in choosing 
the quark type for any computation. We employ the staggered quarks. These have an exact chiral 
symmetry which provides an order parameter for the entire T-^b plane but unfortunately flavour 
and spin symmetry are broken for them on lattice. The existence of the critical point, on the 
other hand, is expected to depend crucially on the number of flavours. Although computationally 
much more expensive. Domain Wall or Overlap Fermions are better in this regard, as they do have 
the correct symmetries for any lattice spacing at zero temperature and density. Introduction of 
chemical potential, ju, for these, however, is not straight-forward due to thek non-locality. Bloch 
and Wettig iQl proposed a way to do this. Unfortunately, it turns out jstl that their prescription 
breaks chiral symmetry . Recently, this problem has been solved [4] and a lattice action with 
nonzero /i and the same chiral symmetries as continuum QCD has been proposed. It will be 
interesting to compare its results with those for the staggered quarks. 

Finite density simulations needed for locating a critical point suffer from another well known 
problem. This one is inherited from the continuum theory itself: the fermion sign problem. 
Assuming Nj flavours of quarks, and denoting by /// the corresponding chemical potentials, the 
QCD partition function is 



Z= Df/exp(-5G) ]~~[DetM(m/,jU/) , 



(1) 



where quarks have been integrated leading to the determinant. The thermal expectation value of 
an observable O is 

/ Df/ exp(-5 c) O Y\f Det M{m /, n /) 
(^) ^ ■ (2) 

Numerical simulations or analytical computations can be done if Det M > for any set of 
{U}. However, Det M is a complex number for any fi 0. This famous phase/sign problem 
is a major stumbling block in extending the lattice techniques to the entire T-fiB plane. Several 
approaches have been proposed in the past decade to deal with it. Let me provide a partial list: 
1) Two parameter Re-weighting [5], 2) Imaginary Chemical Potential |6], 3) Taylor Expansion 
101, 4) Canonical Ensemble method 1^, and 5) Complex Langevin approach jgt]. We employ the 
Taylor expansion approach f?*] to obtain the results discussed in the next section. 



2. Lattice Results 

Our earlier results were obtained by simulating full QCD with two flavours of staggered 
fermions of mass m/Tc = 0.1 on A^, x A??_lattices, with M = 4 and A^, = 8, 10, 12, 16, 24 li^l 
and a finer A^, = 6 with A^,, = 12, 18, 24 | iT]. From the work of the MILC coUaboration 1 12*], we 
know that our lattices correspond to m„/mp - 0.31 + 0.01, leading to a Goldstone pion of 230 
MeV. In order to approach the continuum limit of a — > 0, we proceeded to work on an even finer 
lattice 8 x 32-^, keeping m/Tc =0.1 fixed. As before, the peak of the Polyakov loop susceptibility 
was used to define the critical coupling /3c and plaquette determinations on symmetric lattices 
were used to tune the temperature such that simulation range covered 0.90 < T/Tc < 2.01. Our 
determination of /3c is in excellent agreement with that of Gottlieb et al. (r^. Typically 100- 
200 independent configurations, separated by many autocorrelation lengths were used to make 
measurements of physical quantities. 
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From canonical definitions of number densities «, and susceptibilities Xnu^d' the QCD pres- 
sure P can be seen to have the following expansion in yu: 



AP 

Y4 



P(p, T) P{Q, T) 
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- 2^^""'"" nj\Tl n,\\Tl 



(3) 



where the indices n„ and rid denote the number of derivatives of the partition function with re- 
spect to the corresponding chemical potentials. We set //„ - fid - and m„ = lUd in the 
expressions and construct a series for baryonic susceptibility from this expansion ifioll . Its radius 
of convergence is what we look for 

Successive estimates for the radius of convergence were obtained by using the ratio method 

[r„+i/„+3 = yjnin + \)x^'^^'^ IxT^^T^ ] ^nd the n"' root method [r2/„ = (n\xf lxT''^T"f'"l We 
used terms up to 8th order in yu for doing so. A key point to note is that all coefficients of the 
series must be positive for the critical point to be at real yu, and thus physical. We thus first 
look for this condition to be satisfied and then look for agreement between the two definitions 
above as well as their ^-independence to locate the critical point. The detailed expressions for 
all the terms can be found in where the method to evaluate them is also explained. We 
use stochastic estimators. For terms up to the 8''' order one needs 20 inversions of M on ~1000 
vectors for a single measurement on a given gauge configuration. Work is in progress to double 
these. Our earlier determination of the critical point on A^, = 6 resulted from the constancy for 
both the ratios defined above at T jTc - 0.94, with all the susceptibilities being positive, leading 
Ol in to the coordinates of the endpoint (E) -the critical point-to be /Tc - 0.94 + 0.01, and 
yuf /r^ = 1.8 + 0.1 for that lattice. 
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Figure 1: Comparison of baryon number susceptibility on = 8 with our earlier results on N, =6 and 4 (Left panel). 
Comparison of the expansion coefficients on A', = 6 and 8 lattices. The diagonal solid line displays the ideal case with 
no cut-off elfects. (Right). 

The left panel of Fig. [T] shows a comparison of our new results for the Nt - ^ lattice for 
baryon number susceptibility with those for N, - 6 and 4. The encouraging agreement between 
N, = 2t and 6 suggests that the dimensionless ratios which we employ in all our critical point 
determinations possess only a mild cut-off dependence. The right panel extends this comparison 
to all the susceptibilities we determined for the two lattices with A^, = 8 and 6. The diagonal line 
indicates the trajectory of their expected locations in the ideal case of no lattice cut-off effects. 
This is seen to be so within errors, which still need to be reduced further. 
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Figure 2: The radius of convergence estimates for N, = & (data) along with the N, = 6 critical chemical potential (band) 
at the critical temperature of Te/T^ = 0.94 (Left panel). QCD phase diagram with all known lattice determinations for 
critical point (Right). Our new N, = 8 results, presented in this talk, are denoted by the asterisk. N, = A H fioil and 6 
fnh results are shown by the circles and square respectively. 



Fig. |2]in its left panel displays our new estimates for the A', - 8, obtained by using the two 
methods mentioned above and the coefficients in Fig. [T] The solid band indicates our critical 
chemical potential estimate on A^, = 6 and the temperature chosen is the same as the correspond- 
ing critical temperature Te- We see such behaviour in a small band near this temperature, leading 
to a larger error band on Te, as exhibited in the QCD phase diagram in the right panel along with 
our old results for N, = 6 111 111 and 4 [.10.1 . and those from Budapest- Wuppertal group both [52 of 
which use A^, = 4. 



3. Summary 

The elusive QCD phase diagram in T-fis plane has begun to emerge using first principles 
lattice approach. Our lattice results for A^, = 8 are in very good agreement with those for A^, - 6, 
suggesting the continuum limit to be in sight and the critical point estimate to be robust. 

This work was done on the Blue Gene P of Indian Lattice Gauge Theory Initiative, Tata 
Institute (TIFR), Mumbai. We gratefully acknowledge financial and technical support of TIFR. 
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